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We study scalar, electromagnetic and gravitational test fields in the Hayward, Bardeen and Ayon- 
Beato-Garci'a regular black hole spacetimes and demonstrate that the test fields are stable in all 
these spacetimes. Using the sixth order WKB approximation of the linear ’’axial” perturbative 
scheme, we determine dependence of the quasinormal mode (QNM) frequencies on the characteristic 
parameters of the test fields and the spacetime charge parameters of the regular black holes. We 
give also the greybody factors, namely the transmission and reflection coefficients of scattered scalar, 
electromagnetic and gravitational waves. We show that damping of the QNMs in regular black hole 
spacetimes is suppressed in comparison to the case of Schwarzschild black holes, and increasing 
charge parameter of the regular black holes increases reflection and decreases transmission factor of 
incident waves for each of the test fields. 


I. INTRODUCTION 

It is well known that in reality a black hole (BH) 
can never be in an isolated state. It always interacts 
with matter and fields around it and as the result of 
these interactions it takes a perturbed state. Behaviour 
of the perturbative test scalar, electromagnetic, gravita¬ 
tional or Dirac fields is an interesting topic of the black 
hole physics, of special interest is stability of these fields 
treated as perturbations in linear regime. 

Perturbed BH are characterized by gravitational (elec¬ 
tromagnetic, scalar) waves that are characterized by com¬ 
plex frequency which is called quasinormal mode (QNM). 
Real and imaginary parts of this complex frequency gov¬ 
ern the real oscillation frequency and dissipation or en¬ 
hancing of these oscillations, respectively. If we are able 
to detect these waves we can obtain more precise imag¬ 
inations about BHs. This is why the perturbations of 
a BH are the most momentous topic for taking infor¬ 
mation from the BH. For these purposes, perturbations 
of BHs have been studied for a long time. Metric per¬ 
turbations of the Schwarzschild spacetime were studied 
long time ago by Regge and Wheeler for ’’axial” (or odd 
parity) [l[ and by Zerilli for ’’polar” (or even parity) [2] 
perturbations by adding small perturbation terms onto 
the unperturbed background, requiring that these per¬ 
turbations do not change the stress energy tensor of the 
black hole spacetime. Later Vishveshwara studied nu¬ 
merically scattering of waves on the Schwarzschild BH 
background [3} see for the review Q. Afterwards, Chan¬ 
drasekhar presented monograph about perturbation the¬ 
ory of BHs summarizing all results that had been got 
by that time Jij. All the authors showed that the above 


’Electronic address: b.a.toshmatov@gmail.com 
^Electronic address: ahmadjon@astrin.uz 
t Electronic address: zdenek.stuchlik@fpf.slu.cz 
§ Electronic address: ahmedov@astrin.uz 


mentioned perturbative studies of test fields on station¬ 
ary backgrounds lead to a Schrodinger-like equation with 
complex frequency and a specific effective potential. This 
equation can be solved using semiclassical or numerical 
methods depending on the character of the effective po¬ 
tential. So far, the Schrodinger equation has been solved 
for several potentials in terms of specific functions or nu¬ 
merically [5j], We should note that reviews of QNMs of 
BHs and neutron stars have been written in [6|i8|. Ko- 
dama and Ishibashi extended the calculations of QNMs 
into BHs in higher dimensions in their series of papers 

El- 

One of the standard methods for the QNM calculations 
is the WKB (Wentzel-Kramers-Brillouin) approximative 
method that was applied for the first time by Schutz and 
Will [13. Afterwards, in order to increase accuracy of 
this method it was extended to the third and sixth order 
by Iyer and Will [III and Konoplya [T3| . 

Up to now, QNMs of several BH spacetimes have been 
calculated by several methods, such as the WKB method, 
integration of the wavelike equations, fit and interpola¬ 
tion approaches, Frobenius method, method of contin¬ 
ued fractions, the Mashhoon method and so on fl5l - l4d] ]. 
In [EUIfil ] instabilities of the BHs against to the pertur¬ 
bations have been studied. QNMs of regular and dirty 
BHs has been investigated in 1171 l.j 11 . 

So far, curvature singularity has been considered to be 
the fundamental problem of General Relativity, as the 
curvature singularity cannot be explained by the General 
Relativity itself. As the result of efforts to achieve a new 
solution without singularity, in 1968 Bardeen presented a 
new regular BH solution of modified Einstein field equa¬ 
tions coupled to a nonlinear electrodynamics of a self- 
gravitating magnetic monopole [13 ■ Afterwards, in 1998 
Ayon-Beato and Garcia presented a new charged regu¬ 
lar BH solution based on the Einstein gravity and non¬ 
linear electrodynamics. For finding this regular BH solu¬ 
tion, as a source of charge they took the nonlinear elec¬ 
tric field for the solution of Maxwell field equations [53}. 
All the subtleties of the non-linear electrodynamics and 
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the presence of the regular BH solutions were discussed 
in the series of papers of Bronnikov [53456}. Regular 
BH solutions are discussed also in the framework of the 
so called Modified Gravity (MOG) introduced as tensor- 
vector-scalar (TeVeS) gravitation theory by Moffat [57} . 
In 2006 Hayward found new regular BH solution that is 
free of charge term and quiet similar to the Bardeen one 
in terms of physical aspects [58} . Geodesic motion in the 
field of regular BHs has been studied in [HoHoTt . For the 
regular no-horizon highly curved spacetimes, being com¬ 
plementary to the regular BH spacetimes, the geodesic 
motion has been studied in [62|, the extraordinary ef¬ 
fect of ghost images that do not occur neither in BH 
or naked singularity spacetimes has been exposed in j63j. 
Geodesic structure of the no-horizon spacetimes is of sim¬ 
ilar character as those of the naked singularity Reissner- 
Nordstrom [63467} or Kehagias-Sfetsos spacetimes [ 68 l — 
but differs significantly for the Kerr naked singular¬ 
ity spacetimes I7li 1731 . Here we concentrate on the be¬ 
haviour of the scalar, electromagnetic and gravitational 
perturbative fields in the regular Hayward, Bardeen and 
Ayon-Beato-Garcia (ABG) BH spacetimes. The dynam¬ 
ical stability of black hole solutions in self-gravitating 
nonlinear electrodynamics with respect to arbitrary lin¬ 
ear fluctuations of the metric and electromagnetic field 
has been analyzed in [74} . We calculate the quasinormal 
modes and reflection or transmission coefficients using 
the 6 th order WKB method. 

This paper is organized as follows. Section II gives brief 
description of the Hayward, Bardeen and ABG regular 
BH spacetimes. In section III the perturbative equations 
of the test scalar, electromagnetic and gravitational fields 
in given backgrounds are reduced to the Schrodinger-like 
wave equation. Stability of the perturbative scalar, elec¬ 
tromagnetic and gravitational fields in the regular BH 
spacetimes is treated in the linear regime in section IV. 
Section V is devoted to the calculations of the eikonal 
limit for the frequency of the QNM in the regular BH 
spacetimes under consideration and numerical calcula¬ 
tions are made by the 6 th order WKB approximation. 
In section VI we analyze transmission and reflection of 
incident waves through the potential barrier given by the 
regular BHs. In section VII we summarize our main re¬ 
sults. 


II. REGULAR BLACK HOLES 


Spacetime line element of the spherically symmetric 
regular BH spacetimes can be expressed in the form 


where index i = 1,2,3 corresponds to the Hayward, 
Bardeen and ABG BHs, respectively. 

For the Hayward BH [58| 


TOi(r) = 

for the Bardeen BH fH3l 


Mr 3 


m 2 (r) = 


Mr 3 

(r 2 + Q 2 ) 3/2’ 


for the ABG BH 


( 3 ) 


( 4 ) 


„ , n Mr3 Ql r3 ^ 

3(r) (r 2 +Q 2 ) 3 / 2 2(r 2 + Q 2 ) 2 ' (5) 

Here M, Q 2 and Q 3 are mass, magnetic and electric (or 
magnetic) charge, respectively, while Q± is some real pos¬ 
itive constant. 

The Hayward and Bardeen regular BHs solutions 
asymptotically behave as the Schwarzschild one at large 
distances r —> 00 


lim 

1 - >-00 


/i, 2 (r) = l- — +0 

r 



( 6 ) 


The ABG regular BH solution asymptotically behaves as 
the Reissner-Nordstrom one 


lim / 3 (r) = 1 - — + % + O 

r—>oo T T 



( 7 ) 


When r —> 0, the Bardeen, Hayward and ABG BH space- 
times behave as the de Sitter spacetime 

lim /i, 2 (r) = 1 - + O (r 5 ) , ( 8 ) 

r ~^5 1,2 

9M r 2 

Jim/aM = 1 - + q| + ° ( r ) • (») 

The ABG BH spacetimes may demonstrate an anti de- 
Sitter behaviour near the origin of coordinates, however 
it is related to the no-horizon case (Q > 2 M) |62|. 

For convenience we turn into dimensionless quantities: 
r -A r/M, g -A Q±/M, q -A Q 2 /M and d -A Q 3 /M. 
Then for the Hayward BHs 


to i(r) 



for the Bardeen BHs 


( 10 ) 


ds 2 


-f{r)dt 2 + 


dr 2 + r 2 d0 2 + r 2 sin 2 6dip 2 

/ 0 ) 


(1) 


m 2 (r) 


( r 2 + g 2)3/2 > 


( 11 ) 


with the lapse function f(r) determined by the formula 


and for the ABG BHs 


2 rriiir) 


( 2 ) 


7 *3 

m 3 (r ) = 


fi(r) = 1 


r 


(r 2 + d 2 ) 3 / 2 2 (r 2 + d 2 ) 2 


( 12 ) 
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The horizons of the static spherically symmetric BHs 
are defined by vanishing of the time metric component, 
gtt = 0, i.e., the lapse function f{r 7 qi) = 0. By solving 
this equation with respect to the radial coordinate r, one 
obtains the outer event horizon r + , and the inner hori¬ 
zon r_. One can see from the Fig. |T] that the regular 
Hayward, Bardeen and ABG spacetimes can have two 
horizons, corresponding to the BH spacetimes, one hori¬ 
zon, corresponding to the extremal BH spacetimes, and 
no horizon, corresponding to the no-horizon spacetimes, 
in dependence on the values of the mass m and param¬ 
eters qi. In order to find the critical value of the charge 
q c corresponding to the extremal BH spacetimes (when 
= r_), one can use the system of equations 

gtt = 0, 9tt,r = 0. (13) 

By solving above equations simultaneously, we obtain the 
critical charge of the Hayward spacetimes to be g c ~ 
1.0583 while related radius of the horizon reads r c = 
1.3333. It is shown in Fig. [T| that in the case q < q c 
there are two horizons, in the case q > q c there are no 
horizons. For the Bardeen and ABG spacetimes we arrive 
at the critical charges and related horizon radii expressed 
in dimensionless form: q c ss 0.7698, r c = 1.0887 and 
d c ~ 0.6342, r c = 1.005, respectively. 
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FIG. 1: Radial dependence of lapse function fi(r) of the reg¬ 
ular BHs for different values of the parameter qi. Where qi 
is g, q, d for the Hayward, Bardeen and ABG spacetimes, re¬ 
spectively. 



III. EQUATIONS GOVERNING TEST FIELDS 

Let us consider a massless scalar field $ in the spheri¬ 
cal symmetric regular BH spacetimes, obeying the Klein- 
Gordon equation [zU 

» ee -^= 0 ,, (v=s<ra,«i>) = o, (i4) 

where d L is the partial derivative, and —g is the absolute 
value of the determinant of the spacetime metric ED- 


In order to reduce the Klcin-Gordon equation (fl4l) into 
the two dimensional wave equation, we assume the scalar 
field separated in the standard form 

cj>(t,r,M) = ^M)}T(M) , (15) 

where Yf(0, (f>) is the so called spherical harmonic function 
of degree £ related to the angular coordinates 9, <f>, and it 
fulfills the relations 




1 


■dg (sin Odg) 


1 


sin0 

= -£(£+i)Yr(e,<t>) . 


sin 2 9 


dL Yr(8,0) 

(16) 


Inserting the separated scalar function (11511 into (1141) , we 
obtain the Regge-Wheeler wave equation [lj relating the 
time and radial dependence of the scalar function 


dt 2 


drl 



R(r , t) = 0. 


(17) 


In the Regge-Wheeler equation, r* is the so called ’’tor¬ 
toise” radial coordinate defined by 


r * 



(18) 


Since f(r) is the lapse function governing the event hori¬ 
zon location, r* approaches -oo as r approaches the 
event horizon of the BH from infinity, thus 


r» € (—oo,+oo) ,for r € (r + ,+oo). (19) 


Therefore, the Regge-Wheeler wave-like equation E3 
will be restricted here only to the regions located out¬ 
side the event horizon, r > r + . 

The effective potential V (r) in the expression lH71) de¬ 
pends on the specific field under consideration. For the 
massless scalar fields it reads 


V s (r) 


f(r) 


'£{£ + 1 ) 

r 2 


1 df(r) ' 
r dr 


( 20 ) 


where parameter £ related to the spherical harmonics rep¬ 
resents the orbital angular momentum and takes only 
nonnegative integers. Generalized form of the effective 
potential to higher spin (bosonic) fields can be written in 
the form [H, [501 


V(r) = f(r) 


£(£- 


!) , 2N 2m ( r ) 


+ (l-* 2 ) : 


f 1 df(r) 

2 m(r) ^ 

\r dr 

r 3 J 


( 21 ) 


where the multipole number £ is restricted by £ > s, and 
s is the spin of the perturbative field 


{ 0 , scalar perturbation , 

1 , electromagnetic perturbation , (22) 

2 , gravitational perturbation. 
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FIG. 2: Dependence of V on radial coordinate r for several 
values of orbital quantum number t. 




r/M 


One can deduce from Fig.[2]that the angular momentum 
parameter £ increases the height of the potential barrier 
governed by the effective potential. 

In order to solve the Regge-Wheeler differential equa¬ 
tion (ED we separate the r, t variables in the wave func¬ 
tion assuming 


R(r,t) =ip(r)£(t) . (23) 

Putting the separated function (1731) into (ED, using the 
stationarity of the spacetime, and considering a wave in¬ 
coming from infinity towards the BH, 

m ~ e-* 1 *, (24) 

we obtain the radial part of the Regge-Wheeler wave 
equation in the form of the standard Schrodinger equa¬ 
tion 


+ [w 2 - V{r)] ip{r) = 0 . (25) 

As mentioned before, the QNM frequency can be gener¬ 
ally expressed in the form 

oo = Re(co) + Hm(co) = cor + icoi (26) 

where cor is the real oscillation frequency, and u>i is the 
imaginary oscillation frequency that represents damping 
or enhancing of the oscillatory mode. 


IV. STABILITY OF REGULAR BH 
SPACETIMES 

Equation (1751) is treated considering that there are only 
outgoing waves at infinity and only ingoing waves at the 
event horizon of regular BHs. 

Let us rewrite the wavelike equation (1251) in the form 

- ^ + V(r)ip(.r) - uo 2 ip(r) = 0 . (27) 


In order to study stability of the regular BHs against 
’’axial” perturbations (or, more precisely, stability of the 
perturbative fields in the given spacetimes), we multiply 
both sides of equation (1271) by complex conjugate of the 
wave function ip{r) 

— + V(r)ip(r)ip*(r) — w 2 i/>(7')i/>*(r) = 0 . 

art 

(28) 

Integrating equation (l28l) along whole the range of the 
’’tortoise” coordinate r*, we arrive to the formula 


f* + 00 


V(r) IV’WI 2 + l^'MI 2 dr* - |t~ 

r+oo 

-CO 2 / \ijj{r)\ 2 dr* = 0 ,(29) 


where stands for derivative with respect to the ’’tor¬ 
toise” coordinate r*. 

As mentioned in Sec. Ill, the effective potential related 
to the Hayward, Bardeen and ABG regular BHs, V(r), 
is always positive in the region r* £ (—oo, +oo), i.e., out¬ 
side the event horizon r £ [r + , +oo). Now we rewrite the 
equation (17^1) for the imaginary part using the bound¬ 
ary conditions (EHl) . (1371) and positivity of the effective 
potential outside the horizon [751 ] 


ojr 


|^(+oo)| 2 + |^>(-oo)| 2 + 2 co i 



= 0 . 
(30) 


According to this relation we can deduce that positivity 
of the real part of the frequency, cor > 0, requires the 
imaginary part of the frequency to be negative, coi < 0. 
This means that the perturbative test fields have no expo¬ 
nentially growing modes and the regular BH spacetimes 
are stable under linear odd parity perturbations. 


V. WKB METHOD AND NUMERICAL 
RESULTS 

Boundary conditions. Considering the boundary 
condition for an incoming wave at the event horizon, one 
can see from Fig.[7]that for r —> r+, the effective potential 
V e ff —> 0, therefore, in this limit, the incoming wave 
function ip{r) behaves as 

ip(r) ~ e _ “ r * as r —> r+ (r* -A —oo). (31) 

Considering the boundary condition for an outgoing wave 
at infinity we have to use the fact that all the consid¬ 
ered regular BH spacetimes are asymptotically flat, i.e., 
at infinity they tend to the Minkowski, flat spacetime. 
In the Minkowski spacetime, the effective potential van¬ 
ishes, V e ff = 0. Therefore one can write the boundary 
condition for the outgoing wave as 

ijo{r) ~ e lur * as r — > oo (r* —> + 00 ). (32) 
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WKB method. So far, the Schrodinger equation has 
been solved for several (solvable) potentials in terms of 
special functions or numerically. The general approach to 
solve the Schrodinger equation for properly defined effec¬ 
tive potentials is to reduce this equation to the equation 
for the hypergeometric functions or some other special 
functions. 

One of the methods of solving the Schrodinger equa¬ 
tion is the WKB one. This method was applied for the 
QNMs of test field in the field of BHs for the first time 
by Schutz [12], The WKB method is a semi-analytic 
technique used to solve Schrodinger-type equations such 
as (1251) . It is used mostly for the time-independent case, 
in other words, for an eigenstate of QNM having two 
turning points in the effective potential. However, this 
method also have its validity region. WKB has very good 
accuracy, if the two turning points (solutions of equa¬ 
tion uj 2 — V (r) = 0) are very close to each other. In other 
words, in the case when [V(r) — w 2 ] max -C [w 2 — V(±oo)]. 
In the WKB method, the total energy, w 2 — V(r), is ex¬ 
panded to the Taylor series near the maximum of the 
effective potential nearby the turning points. The WKB 
method has been extended to the third and the sixth or¬ 
der by Iyer and Will 0 and Konoplya 0], respectively. 

In this paper we use the sixth order WKB method of 
calculation of the quasi-normal modes 0 that is gov¬ 
erned by the relation 


i{io 2 -V{r 0 )) 

y/— 2 K"(ro) 


+ X! A i — n + 2 


(33) 


where V"(ro) is the value of the second derivative of the 
effective potential with respect to r at its maximum point 
ro defined by the solution of the equation d,V/dr* | r »=r 0 = 
0 . Aj are constants coming from the second up to the 
sixth order WKB corrections 00 . 

The eikonal limit. It has been already shown 
in [§j, [§] that the WKB method works with large ac¬ 
curacy for large values of the multipole quantum number 
£. For £ » 1, frequency of QNMs of a massless ’’ax¬ 
ial” perturbation field with arbitrary spin s can be found 
analytically by using the first order WKB approach. In 
the case of the Hayward, Bardeen and ABG regular BH 
spacetimes we arrive to the formulae 


ujh = 


I g 3 + r o( r o - 2 ) 

+ r o) 


UJB = 


V + r 2 ) 3 / 2 ~ 2 r 2 
r o(l 2 + r o) 3/2 


iC H (n+ i) 



iC B (n + i) 



(34) 


(35) 


UABG = 


\ 


id 2 + rg ) 2 + r-g [d 2 - 2 yjd 2 + r 2 ^j 


r 2 ( d 2 +r 2 Y 
' + - - iC A BG{n+ -) 


+ °[j) ^ ( 36 ) 


Interestingly, the zeroth order term of the quasinormal 
frequency w, expressed in powers of l/£, does not depend 
on the spin s. Therefore, in any regular BH spacetime 
oscillations of the quasinormal modes with large values 
of the multipole quantum number l have the same fre¬ 
quency for each ’’axial” perturbative test field, indepen¬ 
dent of the spin of the field. 

The cumbersome constants C B , Ch and Cabg are 
functions of the dimensionless parameters g , magnetic 
charge q and electric charge d, respectively. In the vac¬ 
uum cases (g = 0 , q = 0 , d = 0 ), there is Ch = 1 , 
Cb = 1 , Cabg = 1 ■ and the above given expressions are 
identical with the one corresponding to the Schwarzschild 
BH 0. 

For the extremal Hayward BH spacetimes, the effective 
potential reaches its maximum at tq ~ 2.6524 and the 
QNM frequency takes the form 


lim u>h 

(.—^OO 


0.2034 


^+i-* 0 . 8011 (n+i) 


(37) 


In the extremal Bardeen BH spacetimes, there is ro « 
2.3012 and the QNM frequency takes the form 


lim lob ~ 0.2210 


1 

2 


i0.7048(n + —) 


(38) 


In the extremal ABG BH spacetimes, there is r 0 ps 2.1372 
and the QNM frequency takes the form 


lim to abg 
£—>■00 


0.2339 


£+^ -i0.6854(n+i) 


(39) 


The Schrodinger-like wave equation C5l) with the ef¬ 
fective potential (1211) containing the lapse function f(r) 
related to the regular BHs is not solvable analytically. 
Using the sixth order WKB method, we calculate nu¬ 
merically the frequency of QNMs in regular BHs for 
the scalar, electromagnetic and gravitational perturba¬ 
tive fields. It is known that the WKB method demon¬ 
strates high accuracy for low overtones with small imag¬ 
inary part of the QNM frequency if £ > n 0 . 

In Figs [3] - [5] the real and imaginary parts of the QNM 
frequencies are presented for the Hayward, Bardeen and 
ABG regular BHs being given as a function of the WKB 
order. 

In Figs [ 6 ] - [5] dependence of the real and imaginary 
parts of the QNM frequencies are given for the Hayward, 
Bardeen and ABG BHs in dependence on the spacetime 
parameter qi for fixed multipole and overtone l = 2 , n = 0 
parameters. 

It is well known that the l = 0 modes admit only n = 0 
overtones. One can see from the Tab.[3]that for this mode 
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WKB order 



WKB order 


FIG. 3: Real lor (with •) and Imaginary uij (with ■) parts of the QNM frequency as a function of WKB order for scalar 
perturbations (s = 0) of the extremal Hayward (solid), extremal Bardeen (dashed) and extremal ABG (dotdashed) BHs are 
given for the l = 0, n = 0 (left) and l = 3, n = 2 (right) modes. 


(Or, (0[ 



(Or, (Oj 


WKB order 



WKB order 


FIG. 4: Real lor (with •) and Imaginary loi (with ■) parts of the QNM frequency as a function of the WKB order for the 
electromagnetic perturbative fields (s = 1) in the extremal Hayward (solid), extremal Bardeen (dashed) and extremal ABG 
(dotdashed) BHs are given for the l = 1, n = 0 (left) and l = 3, n = 2 (right) modes. 


0J I (Or, OJ, 



WKB order 


FIG. 5: Real lur (with •) and Imaginary un (with ■) parts of the QNM frequency are given as a function of the WKB order 
for the gravitional perturbative fields (s = 2) in the extremal Hayward (solid), extremal Bardeen (dashed) and extremal ABG 
(dotdashed) BHs for the l = 2, n = 0 (left) and l = 3, n = 2 (right) modes. 


the third and the sixth order WKB results demonstrate 
large relative error. For the l = 0, n = 0 mode, such 
an error between the 3th and 6th order WKB approx¬ 
imations has been observed also for the Schwarzschild 

BHs [li¬ 


the QNM frequency related to the Schwarzschild BHs 
(g = 0, q = 0 and d = 0). 


In Fig. [9] intersection of all the curves corresponds to 
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FIG. 6: Dependence of the real (left) ojr and imaginary (right) uji part of the QNM frequencies of the scalar (solid *), 
electromagnetic (dashed •) and gravitational (dotted *) perturbative fields on the charge parameter g of the Hayward BHs is 
given for the l = 2, n = 0 mode. 



FIG. 7: Dependence of the real (left) ojr and imaginary (right) u)i part of the QNM frequencies of the scalar (solid *), 
electromagnetic (dashed •) and gravitational (dotted *) perturbative fields on the magnetic charge q of the Bardeen BHs is 
given for the l = 2, n = 0 mode. 



FIG. 8: Dependence of the real (left) ojr and imaginary (right) u)i part of the QNM frequencies of the scalar (solid *), 
electromagnetic (dashed •) and gravitational (dotted *) perturbative fields on charge d of the ABG BHs is given for the 1 = 2, 
n = 0 mode. 


VI. SCATTERING AND GREYBODY FACTOR 

The Schrodinger-like equation m can be applied to 
scattering of waves in the regular BH spacetimes by us¬ 
ing the standard techniques for tunneling in quantum 
mechanics. The asymptotic solutions of the Schrodinger 
equation ED read: 

ip = A{oj)e~ iuir - + 5(w)e iajr *, r* -► -oo, 

ip = C(w)e- i “ r * + D(to)e iuir *, rw + oo. (40) 

For waves incoming towards the regular BHs from infinity 
there is necessarily B(u>) = 0. The reflection amplitude 


is given by R(u>) = D(ui)/C(uj), while the transmission 
amplitude reads T(w) = A(w)/C(w). Therefore, we can 
write 

ip = T(w)e _I “ r *, r* — oo , 
xP = + R(u)e iur *, r* -a +oo . (41) 

The square of the amplitude of the wave function at a 
particular point of spacetime determines probability of 
finding it in the given point. The wave incoming towards 
a regular BH is partially transmitted and partially re¬ 
flected by the potential barrier. Probability of finding 
this wave in the whole region above the horizon is always 
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FIG. 9: Plots showing the behavior of l = 2, n = 0 QNM frequencies in the complex lj plane for the scalar, electromagnetic 
and gravitational perturbative fields (from left to right) in the Hayward (solid *), Bardeen (dashed ■) and ABG (dotted A) 
BH spacetimes. 


equal to one, therefore the condition 

\R(co)\ 2 + \T(co)\ 2 = l, (42) 

has to be satisfied. Here we should consider following 
three cases: 

(i) oj 2 is much less then the maximum of the effective 
potential V(r 0 ) (to 2 CV(ro)); 

(ii) ui 2 is of the same order as the maximum of the effec¬ 
tive potential V(ro) (w 2 ~ V(r 0 )); 

(iii) oj 2 is much larger then the maximum of the effective 

potential V(ro) (w 2 V(r 0 )). 

In the first case, the greybody factor, i.e., the transmis¬ 
sion coefficient, is close to zero and reflection coefficient 
is nearly equal to one. The third case is opposite to the 
first one, namely the transmission coefficient is close to 
one, while the reflection coefficient is close to zero [79l] . 
The most interesting case is the second one, oj 2 ~ V(r 0 ). 
It is known that the WKB method has high accuracy 
when classical turning points (solution of the equation 
w 2 — V(r 0 ) = 0) are very close to each other. This oc¬ 
curs namely when the condition w 2 ~ V(ro) is satisfied. 
Therefore, in this case we can find the transmission and 
reflection coefficients by applying the sixth order WKB 


approach (see [79] and references therein). Then the re¬ 
flection amplitude reads 


R(u) = (l + e ~ i2nK ^A 


- 1/2 


where 


K(u) 


i (a; 2 - T(r 0 )) 


+ X>. 


2=2 


(43) 


(44) 


From condition (1421) we can find the expression for the 
transmission coefficient that reads 


|2»l =1- 


(l + e ~ l27TK ^ 


- 1/2 


(45) 


In Figs [TO] and |TT] results of numerical calculations of 
the transmission and reflection coefficients for the scalar, 
electromagnetic and gravitational fields with the multi¬ 
pole number l = 2 are shown for the regular Hayward, 
Bardeen and ABG BH spacetimes. 
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VII. CONCLUSIONS 

We have studied QNMs of the linear ’’axial” scalar, 
electromagnetic and gravitational perturbations in the 
Hayward, Bardeen and ABG regular BH spacetimes by 
using the sixth order WKB approximation. Calculations 
have shown that increasing of the spacetime charge pa¬ 
rameter implies monotonic increasing of the real part of 
QNM frequency and monotonic decreasing of the imag¬ 
inary part of QNM frequency, i.e., the damping rate of 
the wave decreases. It means that in the regular BHs os¬ 
cillators are better (slowly damped) than in the field of 
Schwarzschild BHs. In the Hayward and Bardeen regular 
BH spacetimes, damping of the QNMs is always largest 
for the scalar fields, mediate for the electromagnetic fields 
and smallest for the gravitational fields. However, for the 
ABG BH spacetimes the situation is more complex, as 
the scalar fields have damping smaller than the electro¬ 
magnetic fields, if the charge parameter of these space- 
times d > 0.56. 

It has been shown that all of the considered regular 
BHs the linear ’’axial” perturbative scalar, electromag¬ 
netic and gravitational fields are stable. 

The greybody factors of scattering of scalar, elec¬ 
tromagnetic gravitational waves, namely the transmis¬ 
sion and reflection coefficients through the potential bar¬ 
rier (effective potential) of the regular BH spacetimes 
have been calculated numerically under the condition 
uj 2 ~ V(ro) related to the WKB approximation, con¬ 
sidering u> is real. In Figs. fTTTl and fill it is demonstrated 


that for the gravitational perturbative fields probability 
of the wave transmission through the potential barrier is 
larger than for the scalar and electromagnetic ones. In¬ 
crease in the value of the spin of the perturbative fields 
weakens the potential barrier as related to transmission 
of waves through. Furthermore, an increase in charge of 
the regular BHs decreases the transmission and increases 
the reflection of the wave through potential barrier. One 
can conclude that the regular Hayward and ABG BHs 
are favoured in terms of transmission and reflection of 
the wave through their potential barriers, respectively. 
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s=0 

Extremal Hayward (g c = 1.0583) 

Extremal Bardeen ( q c = 0.7698) 

Extremal ABG ( d c = 0.6342) 

overtone, 

multipole 

3rd order WKB 

6tli order WKB 

3rd order WKB 

6tli order WKB 

3rd order WKB 

6tli order WKB 

n=0; 1=0 

0.0894 - 0.0978i 

0.1151 - 0.0768i 

0.1045 - 0.0935* 

0.1211 - 0.0776* 

0.1110 - 0.0958* 

0.1271 - 0.0806* 

n=0; 1=1 

0.3037 - 0.08161 

0.3072 - 0.0828i 

0.3320 - 0.0789* 

0.3343 - 0.0792* 

0.3513 - 0.0812* 

0.3535 - 0.0815* 

n=l; 1=1 

0.2609 - 0.2595i 

0.2716 - 0.2583* 

0.2980 - 0.2468* 

0.3046 - 0.2446* 

0.3181 - 0.2532* 

0.3242 - 0.2511* 

n=l; 1=2 

0.4837 - 0.2486* 

0.4872 - 0.2485* 

0.5337 - 0.2378* 

0.5358 - 0.2375* 

0.5662 - 0.2445* 

0.5680 - 0.2442* 

n=l; 1=3 

0.6950 - 0.2463i 

0.6965 - 0.2463* 

0.7607 - 0.2357* 

0.7615 - 0.2356* 

0.8059 - 0.2424* 

0.8066 - 0.2423* 

n=2; 1=2 

0.4391 - 0.4250* 

0.4432 - 0.4259* 

0.4975 - 0.4042* 

0.4995 - 0.4043* 

0.5307 - 0.4151* 

0.5322 - 0.4153* 

n=2; 1=3 

0.6623 - 0.4159i 

0.6642 - 0.4160* 

0.7342 - 0.3968* 

0.7351 - 0.3968* 

0.7798 - 0.4080* 

0.7805 - 0.4080* 

n=2; 1=4 

0.8769 - 0.4124i 

0.8779 - 0.4124* 

0.9641 - 0.3939* 

0.9646 - 0.3938* 

1.0224 - 0.4050* 

1.0228 - 0.4050* 


TABLE I: Extremal Hayward, Bardeen and ABG BH QNM frequencies given for the scalar perturbative field (s = 0). 


S=1 

Extremal Hayward ( g c = 1.0583) 

Extremal Bardeen ( q c = 0.7698) 

Extremal ABG ( d c = 0.6342) 

overtone, 

multipole 

3rd order WKB 

6tli order WKB 

3rd order WKB 

6tli order WKB 

3rd order WKB 

6tli order WKB 

n=0; 1=1 

0.2630 - 0.0748* 

0.2663 - 0.0784* 

0.2929 - 0.0733* 

0.2957 - 0.0745* 

0.3114 - 0.0758* 

0.3140 - 0.0767* 

CN 

Jl 

O 

II 

s 

0.4851 - 0.0792* 

0.4861 - 0.0795* 

0.5308 - 0.0763* 

0.5314 - 0.0764* 

0.5623 - 0.0787* 

0.5629 - 0.07876* 

n=l; 1=1 

0.2138 - 0.2421* 

0.2277 - 0.2504* 

0.2546 - 0.2316* 

0.2639 - 0.2318* 

0.2740 - 0.2384* 

0.2830 - 0.2377* 

n=l; 1=2 

0.4593 - 0.2417* 

0.4631 - 0.2419* 

0.5103 - 0.2320* 

0.5126 - 0.2318* 

0.5423 - 0.2389* 

0.5443 - 0.2387* 

n=l; 1=3 

0.6780 - 0.2428* 

0.6795 - 0.2428* 

0.7442 - 0.2327* 

0.7451 - 0.2326* 

0.7890 - 0.2395* 

0.7898 - 0.2395* 

n=2; 1=2 

0.4134 - 0.4140* 

0.4183 - 0.4153* 

0.4732 - 0.3947* 

0.4757 - 0.3949* 

0.5058 - 0.4059* 

0.5077 - 0.4062* 

n=2; 1=3 

0.6449 - 0.4101* 

0.6469 - 0.4102* 

0.7175 - 0.3919* 

0.7185 - 0.3919* 

0.7627 - 0.4032* 

0.7635 - 0.4032* 

n=2; 1=4 

0.8635 - 0.4088* 

0.8646 - 0.4088* 

0.9513 - 0.3909* 

0.9517 - 0.3908* 

1.0092 - 0.4021* 

1.0096 - 0.4021* 


TABLE II: Extremal Hayward, Bardeen and ABG BH QNM frequencies are given for the electromagnetic perturbative fields 
(« = !)■ 


s=2 

Extremal Hayward ( g c = 1.0583) 

Extremal Bardeen ( q c = 0.7698) 

Extremal ABG ( d c = 0.6342) 

overtone, 

multipole 

3rd order WKB 

6tli order WKB 

3rd order WKB 

6tli order WKB 

3rd order WKB 

6tli order WKB 

n=0; 1=2 

0.3944 - 0.0735* 

0.3954 - 0.0753* 

0.4321 - 0.0714* 

0.4328 - 0.0719* 

0.4689 - 0.0818* 

0.4698 - 0.0819* 

CO 

Jl 

o 

II 

0.6342 - 0.0779* 

0.6347 - 0.0780* 

0.6916 - 0.0750* 

0.6919 - 0.0751* 

0.7374 - 0.0803* 

0.7377 - 0.0803* 

n=l; 1=2 

0.3620 - 0.2256* 

0.3675 - 0.2349* 

0.4083 - 0.2180* 

0.4100 - 0.2209* 

0.4445 - 0.2505* 

0.4473 - 0.2495* 

n=l; 1=3 

0.6156 - 0.2359* 

0.6174 - 0.2360* 

0.6767 - 0.2267* 

0.6778 - 0.2266* 

0.7222 - 0.2425* 

0.7232 - 0.2424* 

n=l; 1=4 

0.8422 - 0.2394* 

0.8430 - 0.2394* 

0.9210 - 0.2297* 

0.9214 - 0.2296* 

0.9788 - 0.2412* 

0.9792 - 0.2412* 

n=2; 1=2 

0.3037 - 0.3901* 

0.3218 - 0.4190* 

0.3657 - 0.3736* 

0.3656 - 0.3888* 

0.4030 - 0.4290* 

0.4043 - 0.4294* 

n=2; 1=3 

0.5803 - 0.3992* 

0.5829 - 0.3997* 

0.6484 - 0.3823* 

0.6497 - 0.3823* 

0.6935 - 0.4090* 

0.6945 - 0.4091* 

n=2; 1=4 

0.8155 - 0.4025* 

0.8167 - 0.4025* 

0.8994 - 0.3854* 

0.9000 - 0.3854* 

0.9572 - 0.4047* 

0.9577 - 0.4047* 


TABLE III: Extremal Hayward, Bardeen and ABG BH QNM frequencies are given for the gravitational perturbative fields 
(s = 2). 
























































12 


s = 0 

Hayward 

Bardeen 

ABG 

9/9c, q/qc, d/d c 

3rd order WKB 

6th order WKB 

3rd order WKB 

6th order WKB 

3rd order WKB 

6th order WKB 

0.1 

0.2911 - 0.0980* 

0.2929 - 0.0977* 

0.2914 - 0.0979* 

0.2932 - 0.0977* 

0.2915 - 0.0980* 

0.2933 - 0.0977* 

0.2 

0.2912 - 0.0979i 

0.2930 - 0.0977* 

0.2924 - 0.0977* 

0.2941 - 0.0975* 

0.2928 - 0.0979* 

0.2946 - 0.0976* 

0.3 

0.2915 - 0.0977i 

0.2933 - 0.0975* 

0.2939 - 0.0973* 

0.2957 - 0.0971* 

0.2949 - 0.0977* 

0.2967 - 0.0975* 

0.4 

0.2919 - 0.0973i 

0.2938 - 0.0972* 

0.2962 - 0.0966* 

0.2980 - 0.0966* 

0.2980 - 0.0974* 

0.2998 - 0.0972* 

0.5 

0.2927 - 0.0966* 

0.2947 - 0.0966* 

0.2993 - 0.0957* 

0.3011 - 0.0957* 

0.3022 - 0.0969* 

0.3039 - 0.0968* 

0.6 

0.2939 - 0.0955i 

0.2960 - 0.0957* 

0.3033 - 0.0945* 

0.3051 - 0.0946* 

0.3077 - 0.0961* 

0.3094 - 0.0961* 

0.7 

0.2956 - 0.0938i 

0.2979 - 0.0943* 

0.3083 - 0.0927* 

0.3102 - 0.0929* 

0.3148 - 0.0949* 

0.3165 - 0.0950* 

0.8 

0.2978 - 0.0913* 

0.3005 - 0.0920* 

0.3146 - 0.0900* 

0.3166 - 0.0903* 

0.3239 - 0.0928* 

0.3257 - 0.0930* 

0.9 

0.3006 - 0.0875i 

0.3038 - 0.0885* 

0.3225 - 0.0859* 

0.3247 - 0.0862* 

0.3358 - 0.0890* 

0.3378 - 0.0893* 

0.99 

0.3034 - 0.0823i 

0.3069 - 0.0835* 

0.3310 - 0.0798* 

0.3333 - 0.0801* 

0.3496 - 0.0823* 

0.3518 - 0.0826* 

1.0 

0.3037 - 0.0816i 

0.3072 - 0.0828* 

0.3320 - 0.0789* 

0.3343 - 0.0792* 

0.3513 - 0.0812* 

0.3535 - 0.0815* 


TABLE IV: Dependence of the Hayward, Bardeen and ABG BH QNM frequencies on the parameters qi/q c are given for the 
scalar perturbative fields (s = 0). 


|7f |7f 17’1- 











FIG. 10: Transmission coefficients for the scalar, electromagnetic and gravitational (from left to right) fields in the Hayward, 
Bardeen and ABG BH spacetimes (from top to bottom). There is i = 2, q = 0 (solid), q = 0.75g c (dashed), q = 0.9g c 
(dotdashed), q = q c (dotted). 
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FIG. 11: (color online). Reflection coefficients for the scalar, electromagnetic and gravitational (from left to right) fields in the 
Hayward, Bardeen and ABG BH spacetimes (from top to bottom). There is i = 2, q = 0 (solid), q = 0.75 q c (dashed), q = 0 .9q c 
(dotdashed), q = q c (dotted). 




































